基于内聚力模型的快速荷载下CFRP-混凝土界面剪切剥离模拟

鉴于纤维增强复合材料(Fiber Reinforced Polymer, FRP)已在混凝土结构加固中普遍应用, 能有效提高结构的承载和变形能力. 为揭示动载作用下FRP-混凝土界面的剥离机理, 本文基于ABAQUS平台, 采用内聚力模型模拟CFRP-混凝土界面层, 实现了快速荷载下CFRP-混凝土界面剥离的高效模拟. 结果表明: CFRP表面应变在加载过程中由加载端向自由端方向传递, 随着加载速率的提高, 界面承载能力也随之提高; 界面峰值剪应力也存在显著的应变率效应; 模拟结果与试验结果基本吻合, 说明了新方法的有效性.     

梁柱简化模型在自立塔塔线体系动态响应中的适用性研究

对输电塔进行合理简化可以提高塔线体系动力学仿真的效率。本文给出自立塔梁柱简化模型的计算方法，并提出利用梁柱简化模型计算方法建立自立塔塔线体系整体模型，同时采用桁梁混合模型建立精细化塔线体系整体模型，对两种模型塔线体系静力特性及振型和固有频率等动力特性进行对比分析。以脱冰工况为例，采用生死单元技术将施加在输电线节点上的集中质量单元杀死来模拟脱冰，实现对塔线体系动力学响应的有限元模拟，研究塔线体系简化模型在动态响应中的适用性。结果表明，两种模型弯曲变形误差小，低阶的振型相同，固有频率值误差小，动力特性基本相同；脱冰工况下，自立塔节点位移和塔材内力时程曲线一致，在提高计算效率的情况下，能有效保证计算精度。     

Versatile anisotropic mesh adaptation methodology applied to pure quantity of interest error estimator. Steady,laminar incompressible flow

We introduce a new flexible mesh adaptation approach to efficiently compute a quantity of interest by the finite element method. Efficiently, we mean that the method provides an evaluation of that quantity up to a predetermined accuracy at a lower computational cost than other classical methods. The central pillar of the method is our scalar error estimator based on sensitivities of the quantity of interest to the residuals. These sensitivities result from the computation of a continuous adjoint problem. The mesh adaptation strategy can drive anisotropic mesh adaptation from a general scalar error contribution of each element. The full potential of our error estimator is then reached. The proposed method is validated by evaluating the lift, the drag, and the hydraulic losses on a 2D benchmark case: the flow around a cylinder at a Reynolds number of 20.     

A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

A note on 3‐partite graphs without 4‐cycles

Let $C_4$ be a cycle of order 4. Write $ex\left(n,n,n,C_4\right)$ for the maximum number of edges in a balanced 3‐partite graph whose vertex set consists of three parts, each has $n$ vertices that have no subgraph isomorphic to $C_4$. In this paper, we show that $ex\left(n,n,n,C_4\right)\ge \frac{3}{2}n\left(p+1\right)$, where $n=\frac{p\left(p?1\right)}{2}$ and $p$ is a prime number. Note that $ex\left(n,n,n,C_4\right)\le \left(\frac{3\sqrt{2}}{2}+o\left(1\right)\right)n^{\frac{3}{2}}$ from Tait and Timmons's works. Since for every integer $m$, one can find a prime $p$ such that $m\le p\le \left(1+o\left(1\right)\right)m$, we obtain that $\underset{n\to \infty }{\lim }\frac{ex\left(n,n,n,C_4\right)}{\frac{3\sqrt{2}}{2}{n}^{\frac{3}{2}}}=1$.     

安培环路定理的拓扑学描述

近年来拓扑学在量子力学中得到了广泛的运用.本文将安培环路定理积分式重新表达为一矢量场在轮胎参数面上的第一类陈数积分.数值模拟展示了该积分值为一整数即第一陈数,其代表矢量场的整体性质:当经历连续变换时,矢量场的局部数值发生改变但整体积分值即陈数仍保持不变;若陈数发生改变,则表明矢量场变换的连续性条件发生破坏,矢量场出现奇点.进一步通过高斯映射将该矢量场从参数轮胎面映射到单位球面上,并给出了第一陈数的直观几何意义.理论和数值结果揭示了安培环路定理的拓扑学本质,表明拓扑概念在经典物理学中也会有广泛应用.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.